Experimental observation of spatial antibunching of photons 
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We report an interference experiment that shows transverse spatial antibunching of photons. Using 
collinear parametric down-conversion in a Young-type fourth-order interference setup we show in- 
terference patterns that violate classical Schwarz inequality and should not exist at all in a classical 
description. 



Photon antibunching in a stationary field is recognized 
as a signature of nonclassical behavior, for its description 
is not possible in terms of a nonsingular positive Glauber- 
Sudarshan P distribution |ffl]. It is well known that any 
state of the electromagnetic field that has a classical ana- 
log can be described by means of a positive P distribution 
which has the properties of a classical probability func- 
tional over an ensemble of coherent states. 

The classical intensity correlation function for station- 
ary fields must obey the following inequality El: 



(I(r,t)I(v,t + T)) < (I 2 (r,t)). 



(1) 



All field states described in terms of a positive nonsin- 
gular P distribution must obey the standard quantum 
mechanical counterpart of ([!]), where products of inten- 
sities are replaced by ordered products of photon density 
operators jjj, that is, 



(T:i(v,t)i(r,t + r) :) < <:i 2 (r,i) :), 



(2) 



where T : : stands for time and normal ordering. Photon 
density operators are defined as 



J(r,t)=Vt(r,t)V(r,t), 



(3) 



where 



V(r,t)=^a k , CT e k , CT e J ( kr -^ 

k,(r 



o-k.a is the annihilation operator for the mode with wave 
vector k and polarization a, e^ a is the unit polarization 
vector, and u> — ck. 

Inequality (|J) means that for such class of fields, pho- 
tons are detected either bunched or randomly distributed 
in time. Photon antibunching in time, characterized by 
the violation of (|J), was predicted by Carmichael and 
Walls || , Kimble and Mandel || , and was first observed 
by Kimble, Dagenais and Mandel in resonance fluores- 
cence Q. 

Let us now turn to space domain and consider that the 
transverse field profile of a given stationary light beam 
propagating along z direction is described by a complex 
stochastic vector amplitude V(p, t) with an associated 



probability functional V(V). Here, p lies in a plane trans- 
verse to the propagation direction. The average intensity 
at a point p is 

(J(p,t)> = <V*(p,t)V(p,t)> = J V(V)\V(p,t)\ 2 dV, 

(4) 

and the two-point intensity correlation function 

rW( Pl , P2 ,T) = (i( Pl ,t)i( P2 ,t + T)) 

is 

r^( Pl ,p 2 ,T) = J v(v)\v(p u t 1 )\ 2 \v(p2,t 2 )\ 2 dv. 

(5) 

Its time dependence is restricted to the difference r = 
t\ —t 2 , since the field is assumed to be stationary. In the 
space domain, the concept analogous to stationarity is 
homogeneity. For a homogeneous field, the expectation 
value of any quantity that is a function of position is 
invariant under translation of the origin [hi . In particular, 



and 



r( 2 ' 2 )( Pl ,p 2 ,r) = r^ 2 )(^,r) 



(I«(p + 6,t + T)) = (I«(p,t)), 



(6) 
(7) 



where S = p\ — p 2 and N = 1, 2, . 
Applying Schwarz inequality, 



{I(p,t)I(p + 6,t + T)) 2 < (I 2 (p,t))(I 2 (p + 8,t + T)). 

(8) 

By means of ([?]), 

(I(p,t)I(p + 8,t + r)) < (I 2 (p,t)). (9) 
Quantum mechanically, 

(T:i(p,t)i(p + S,t + r) :) < (: I 2 (p,t) :), (10) 
that is, 

r (2 ' 2) (^^) < r (2 ' 2) (o,o). (li) 



Analogously to what was concluded from inequality 
(||), for field states represented by positive nonsingular 
Glauber-Sudarshan distributions, that is, fields that ad- 
mit the classical stochastic description assumed above, 
inequality ([To]) implies that photons are detected either 
spatially bunched or randomly spaced in a transverse de- 
tection screen. Spatial antibunching of photons has been 
predicted by some authors and a possible exper- 

iment was recently proposed to observe it in squeezed 
states |,|. 

In this work we show that strong antibunching in one 
transverse direction can be observed in down-converted 
light, violating ( |l0|) by several standard deviations. The 
effect is produced by fourth-order interference of a two- 
photon beam diffracted by a birefringent double-slit. 
The experimental setup is depicted in Fig. 1. A light 
beam of A = 702 nm is produced by collinear type II 
down-conversion in a 2 mm-long nonlinear crystal (BBO) 
pumped by an Argon laser beam with A = 351 nm. The 
u. v. beam transmitted by the crystal is removed from the 
down-converted beam by a laser mirror (M) transparent 
to 702 nm. A birefringent double slit (S) is constructed 
as follows. A single slit of dimensions 0.60mmx5mm is 
divided in two by a 0.20 mm wide absorbing strip, defin- 
ing two parallel slits of dimensions 0.20mmx5mm. In 
front of each slit there is a quarter wave plate (Qi and 
Q2), as shown in Fig. 2. One wave plate (Qi) has its 
fast axis parallel to the slits, whereas the other one (Q2), 
has its fast axis perpendicular to the slits. With such 
alignment, the waveplates introduce a phase difference 
of 7r between the two slits. This arrangement is placed in 
the down-converted beam, 38 cm from the crystal. The 
pumping beam is focused right on the plane of the double 
slit by a lens (L) of 500 mm focal length. Assuming that 
the beams are propagating along the z direction, this 
focusing causes the fourth-order correlation function of 
the down-converted beam to be concentrated on points 
satisfying £1 + £2 = 0, where £1 and £ 2 are position 
vectors on the plane of the double slit po[ . The focus- 
ing is essential to produce the appropriate spatial de- 
pendence of the diffracted field p| . In order to make 
possible that two detectors (Di and D2) share the same 
transverse position without being limited by their phys- 
ical dimensions, a beam splitter (BS) is inserted in the 
down-converted beam, with D\ and D2 placed in front of 
each exit port. In front of each detector, there is a sin- 
gle slit of dimensions 0.20mmx3mm aligned horizontally 
(parallel to the slits in S) , followed by an interference fil- 
ter with a bandwidth of 40 nm, centered at 690 nm, and 
a lens focused on the detector's active area. The opti- 
cal path length from the double slit S to the detectors 
D\ and D2 is 70 cm. D\ and D2 are mounted on pre- 
cision translation stages and their vertical positions are 
set by computer-controlled stepping motors. Single and 
coincidence counts were measured while D\ and D2 were 
scanned in the vertical direction (x axis), as will be de- 



scribed below. 

Ideally, the transverse fourth-order correlation func- 
tion r( 2,2 )(p, 0) is proportional to the coincidence rate be- 
tween two punctual detectors separated by p, with a neg- 
ligible resolving time. Since the detectors are not punc- 
tual and the coincidence resolving time is finite (10 ns 
in our setup), what was actually measured is a convolu- 
tion of r( 2,2 )(<5, r) with the sampling window Ax Ay At, 
where Aa; and Ay represent the dimensions of the detec- 
tor entrance slit (0.20mmx3mm) and At is the resolving 
time of the coincidence counter (10 ns). For the purpose 
of demonstrating the effect, however, we will ignore this 
correction by considering Ax ~ 0, Ay — > 00, and At ~ 0. 
Under these conditions, it is possible to show |ll[] that for 
small displacements, the coincidence rate is proportional 
to 



cos 



2nd. 
-Xz~ {X2 



Xl) 



(12) 



where A is the wavelength of the down-converted field 
(702 nm), d is the double slit separation (0.40mm), z is 
the optical path length between the double slit and the 
detectors (70 cm), X\ and X2 are the vertical positions of 
detectors D\ and D2, respectively. 

Before taking correlation measurements, the accuracy 
of vertical positioning was checked by the following pro- 
cedure. With the double slit S removed, a horizontally 
aligned wire was stretched in front of the beam splitter, 
at x = 0, and single counts were registered in sampling 
times of 5 s, while the detectors were scanned vertically. 
The result is shown in Fig. 3. The two counting profiles 
are not identical due to differences in the overall quantum 
efficiencies of D\ and L*2- 

Figures 4 to 7 summarize the results of single counts 
and coincidence measurements taken in sampling times 
of 1000 s in several different situations. All coincidence 
patterns were fit to expression ( |l2] ) plus a background. 
Vertical error bars are statistical with two standard de- 
viations in length, whereas horizontal ones correspond to 
the width of detectors entrance slits. The results shown 
in Fig. 4 refer to the situation in which detector D2 is 
kept at X2 — and detector D\ is scanned vertically. The 
single counts, although not constant, do not show any 
oscillation to which one could attribute the oscillation in 
coincidences. The same is true in Fig. 5, where detector 
D\ was kept in x\ = and D2 was scanned vertically. 
When D\ and D2 were scanned together (xi = X2), a 
fairly constant background of coincidences was recorded, 
as shown in Fig. 6. This background, that should be zero 
as well as the minima in Fig. 4 and Fig. 5, is due to the 
finite width of the detectors entrance slits (0.20 mm). A 
final measurement was performed by scanning D\ with 
D2 kept in the position X2 — —0.55 mm, which corre- 
sponds to a maximum of coincidences in Fig. 5. The 
results are plotted in Fig. 7, showing that the minimum 
in coincidences was displaced to x\ — X2 — —0.55 mm. 



All the interference patterns shown here satisfy 

r (2 > 2 \s,o) > r (2 ' 2) (o,o), (13) 



in a clear violation of expression (11), characterizing the 
presence of transverse spatial antibunching of photons. 

Let us analyze these results from another point of view. 
Some years ago, it was pointed out |L2| that all second- 
and fourth-order optical interference effects observed so 
far have close classical analogs with the same harmonic 
pattern, differing only in their visibilities. This is not 
the case of our results. Since the minimum of fourth- 
order interference occurs for X\ = x-i in the absense of 
second-order interference, any classically predicted visi- 
bility different from zero would violate Schwarz inequal- 
ity. Therefore, our results can be regarded as a truly 
quantum fourth-order interference effect. 
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FIG. 1. Experimental setup. L is a lens of focal length 
/ = 500 mm, BBO is a 2mm-long /3-BaB204 nonlinear crys- 
tal cut for collinear type II 351 nm^702 nm down-conversion, 
M is a u. v. high reflectance mirror, Q\ and Q2 are quar- 
ter-wave plates, S is a double slit, BS is a 50:50 beam split- 
ter, Di and D2 are avalanche photo-diodes working in photon 
counting mode. 
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FIG. 2. The birefringent double slit. Qi and Q2 are quar- 
ter-wave plates aligned with orthogonal fast axes, and S is a 
double slit with clear apertures of 0.20mmx5mm separated 
by a 0.20 mm obstacle. 
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FIG. 3. Di(o) and I?2(«) single counts taken with a 
0.20 mm diameter wire stretched horizontally in front of the 
beam splitter and the double slit removed. This measurement 
was taken in order to check the accuracy of detectors vertical 
positioning. 
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FIG. 4. Single counts (<>) and coincidences (•) taken with 
Di kept in X2 = and Di scanned vertically. 
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FIG. 6. Di single counts (o) , D2 single counts (o), and 
coincidences (•) taken when both D\ and D2 were scanned 
vertically, keeping x\ — X2- 
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FIG. 7. Single counts (o) and coincidences (•) taken with 
Z?2 kept in X2 = —0.55 mm and Di scanned vertically. 
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FIG. 5. Single counts (o) and coincidences (•) taken with 
Di kept in x\ = and D2 scanned vertically. 



